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ABSTRACT 

This paper discusses the development 
of a multilayer interference analysis 
and simulation program which is used 
to evaluate interference between non- 
geostationary and geostationary satel- 
lites. In addition to evaluating 
interference, this program can be used 
in the development of sharing criteria 
and coordination among various 
Mobile Satellite Services(MSS)sys- 
tems. A C + + /Windows implementa- 
tion of this program, called Globalstar 
Interference Simulation Program 
(GISP), has been developed and will 
be demonstrated in the conference. 

INTRODUCTION 

Since the relative position and pointing 
of the Geo-Stationary Orbit(GSO) and 
Low Earth Orbit(LEO) satellites con- 
stantly change, the parameters affect- 
ing interference also change with time 
and satellite positioning. With these 
changing parameters, the interference 
situation cannot be evaluated fully 
using conventional techniques of inter- 
ference analysis. To gain better under- 
standing of the interference situation 
and to develop sharing criteria, this 


multilayer interference analysis and 
simulation program was developed. 

In order to help visualize the complex 
interference environment of MSS sys- 
tems, we use the analogy of utilization 
of a single spectrum analyzer to view 
the interference signals. This spectrum 
analyzer can be connected to a ter- 
minal of either a victim satellite or a 
victim user antenna. 

INTERFERENCE TO A SATELLITE 

The total interference signal received 
by a victim satellite antenna terminal 
measured by the fictitious spectrum 
analyzer can be represented by the fol- 
lowing summations. 
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Where 

• is the total power 
spectral density at the ter- 
minal of kth antenna beam of 
victim satellite i, at time t and 
frequency f. It is measured in 
dBW/Hz. 

• N s (t) is the number of inter- 
fering satellites visible to vic- 
tim satellite i at time t. 

• Nbj is the number of multiple 
beams of interfering satellite j. 

• Gijc(Oij(t),<pij(t )) is the gain of 
the kth antenna beam of vic- 
tim satellite i at the angle 
6ij(t) and <pij(t), and 6ij(t) and 
<pij(t) are the corresponding 
victim satellite i’s local spheri- 
cal coordinate angles towards 
interfering satellite j at the 
time instant t. 

• Gjj,(Oji(t),<pji(t )) is the gain of 
the pth antenna beam of inter- 
fering satellite j at the angle 
6ji(t) and <pji(t), and 6ji(t) and 
<Pji(t) are the corresponding 
interfering satellite j’s local 
spherical coordinate angles 
towards victim satellite i at 
the time instant t. 

• Pojp is the power spectral den- 
sity into the terminal of the 
pth antenna beam of interfer- 
ing satellite j. 

• Lij(t) is the path loss between 
victim satellite i and interfer- 
ing satellite j at time t. 


In the second summation, 

• Nu(t ) is the number of inter- 
fering mobile user units 
visible to victim satellite i at 
time t. 

• Gijc(Oij(t),<pij(t)) is the gain of 
the kth antenna beam of vic- 
tim satellite i at the angle 
6ij(t) and <pij(t), and 0ij(t) and 
<pij(t) are the corresponding 
victim satellite i’s local spheri- 
cal coordinate angles towards 
interfering user unit j at the 
time instant t. 

• Gj0ji(t),<pji(t )) is the gain of 
the antenna beam of interfer- 
ing mobile user j at the angle 
Qji(t) and <Pji(t), and Oji(t) and 
<pji(t) are the corresponding 
interfering user j’s local 
spherical coordinate angles 
towards victim satellite i at 
the time instant t. 

• Poj is the power spectral den- 
sity into the terminal of the 
antenna beam of interfering 
mobil unit j. 

• Lij(t) is the path loss between 
victim satellite i and interfer- 
ing mobile user unit j at time t. 


In the above formula, the first summa- 
tion corresponds to the interference 
from all other satellites, whereas the 
second summation corresponds to inter- 
ference from all ground based mobile 
user units. 
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INTERFERENCE TO A MOBILE 
UNIT 

There are essentially two types of inter- 
ference to a mobile user unit; the inter- 
ference from all satellites who transmit 
in the same band as the user receiver, 
and other mobile users. This program 
currently simulates the first kind and 
the future plans include the implemen- 
tation of second kind. To measure the 
total satellite interference to a mobile 
user unit, the fictitious spectrum 
analyzer is connected to terminal of an- 
tenna of the victim mobile user unit. 
The interference that the spectrum 
analyzer measures can be represented 
by the following formula. 
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Where 

• Pj(tf) is the total power 
spectral density at the ter- 
minal of antenna of victim 
user j at time t and frequency 
f. It is measured in dBW/Hz. 

• N s (t) is the number of inter- 
fering satellites visible to vic- 
tim mobile user unit j at time 
t. 

• Nbi k the number of multiple 
beams of interfering satellite i. 

• Gijc(6ij(t),<pij(t )) is the gain of 
the kth antenna beam of inter- 
fering satellite i at the angle 
0y(O and <pij(t), and %(/) and 


<pij(t) are the corresponding 
interfering satellites i’s local 
spherical coordinate angles 
towards victim mobile user j 
at the time instant t. 

• Gj(Oji(t),<pji(t )) is the gain of 
antenna beam of victim 
mobile user j at the angle 
6ji(t) and <Pji(t), and 6>y,(f) and 
<Pji(t) are the corresponding 
victim mobile user j’s local 
spherical coordinate angles 
towards interfering satellite i 
at the time instant t. 

• Poik is the power spectral den- 
sity into the terminal of the 
kth antenna beam of interfer- 
ing satellite i. 

• Lij(t) is the path loss between 
victim user unit j and interfer- 
ing satellite i at time t. 

MULTILAYER MODELING OF IN- 
TERFERENCE 

A systematic approach to partition the 
complex interference situation among 
multiple MSS satellite systems (both 
GSO and non-GSO) is to develop a 
multilayer interference model such that 
a computer analysis and simulation pro- 
gram can be developed around each 
layer. The formulas just presented are 
at the heart of GISP. In the following 
sections, major layers of the program 
are discussed. 

GEOMETRIC LAYER 

Since the relative position and pointing 
of GSO and LEO satellites constantly 
change, it is necessary to develop some 
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common reference systems for all satel- 
lite systems such that all satellites and 
user terminals can have a specific coor- 
dinate with respect to the same refer- 
ence point. In this layer, the center of 
the earth has been identified as the 
center of this reference system as 
shown in Figure 1. Each interference 
source (i.e., a transmitter) and victim 
(i.e., a receiver) is identified by its own 
coordinate with respect to the center 
of the earth. The z-axis of the local 
coordinate system of a satellite or a 
user positioned at any point is defined 
as the vector from this particular point 
to the center of the earth. The y axis is 
defined in the negative direction of the 
normal to the plane of the orbit. A 
right hand coordinate system dictates 
that the x axis be in the plane of the 
satellite velocity . 

The main purpose of this geometric 
layer is to determine whether an inter- 
ference source is visible by the victim 
or otherwise shadowed by earth. Fur- 
ther processing is carried out to deter- 
mine the level of interference if and 
only if the victim is in the line of sight 
of the interferer. Scattering reflection 
and refraction effects are considered 
secondary effects and are not included 
in the program at this stage. Major 
parameters to be input into this 
geometric layer include :satellite al- 
titude, inclination angle, orbit phasing, 
orbit location or sub-satellite point (for 
GSO), orbit eccentricity, apogee and 
perigee, and others. 

The computer program simulates or- 
bits and constellations of different 
GSO or LEO systems, calculates the 
relative positions of each point (either 
an interference source or victim) and 
determines whether two points (i.e. 


source and victim) are in line-of-sight 
or are blocked by the earth at any time 
of the simulation period. 

ANTENNA LAYER 

Once the line-of-sight between the in- 
terferer and the victim is established, 
the interaction between the source an- 
tenna and the victim antenna has to be 
modeled. Figure 2 illustrates the inter- 
action among a single interference 
source and multiple victims. At this 
point, the local coordinate system of 
the satellite or user is defined and the 
pattern angle s 6 and <p are calculated 
relative to these local coordinate sys- 
tems. These angles are then used to 
determine the antenna gain levels of 
interferer toward the victim and the an- 
tenna gain of victim toward the interfer- 
ing source. The relative distance and 
spatial loss be tween the source and vic- 
tim is calculated by the geometric layer 
processor of the program. Thus the 
relative gain between the source and 
victim along the line-of-sight direction 
can be determined. 

MOBILE USER UNITS 

For simplicity, a uniform model of user 
distribution is used in this program. 

The earth is divided into cell regions 
bounded by certain latitudes and lon- 
gitudes. Several mobile user units are 
then lumped at the center of each cell. 
This includes users of various MSS sys- 
tems. In the uniform distribution 
model, it is assumed that users of dif- 
ferent MSS systems are distributed 
evenly and uniformly over the same 
geographic area. This situation is il- 
lustrated in Figure 3. 
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By combining the traffic distribution 
profile of a MSS satellite system and 
the user distribution, interference be- 
tween users and satellites, and between 
satellite and satellite can be simulated 
for a long period of time. To reduce the 
complexity of the problem, user-to- 
user interference is not currently con- 
sidered. 

PROGRAM OUTPUT 

The program "measures" the total inter- 
ference Power Spectral Density (PSD) 
at the terminal of the victim’s antenna. 
The scope of the spectram analyzer is 
displayed graphically by the prgram. As 
the program steps in time, the location 
of the victim is also displayed graphical- 
ly on a flat earth. The interference 
power information and orbital 
dynamics information are captured into 
separate files for further analysis. This 
information can be examined directly 
from the program or from any text 
editor. 


To make the program output more 
comprehensive, several data reduction 
and statistical programs were 
developed. The results of the statistical 
analysis is also captured into a file that 
can be accessed from the program. The 
statistical and power spectral density in- 
formation can be used to estimate the 
aggregate Io and its impact on opera- 
tion of various MSS systems. 

Eb/(No + Io) for any scenarios and any 
MSS system can be calculated to 
evaluate the system operation and link 
integrity and circuit capacity. Further- 
more, the information provided by 
GISP can be used to allocate expected 
interference noise from different MSS 
systems, and thus establish a criteria 
for frequency sharing among various 
MSS systems. 

GISP has been developed through the combined 
efforts of team members Ernie Aylaian, Chuck 
GantZy Steve Gileno, and Murray Steinberg . Their 
contributions are gratefully acknowledged . 



Figure 1. Illustration of various coordinate systems 
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Figure 2. Interference among various satellites 



Figure 3. Illustration of user distribution cells 
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